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Abstract 

First-principles thermodynamic models based on the cluster expansion formalism, monte-carlo sim- 
ulations and quantum-mechanical total energy calculations are employed to compute short-range-order 
parameters and diffuse-antiphase-boundary energies in hep-based a- Ti-Al alloys. Our calculations un- 
ambiguously reveal a substantial amount of SRO is present in a-Ti-6 Al and that, at typical processing 
temperatures concentrations, the DAPB energies associated with a single dislocation slip can reach 25 
mj/m 2 . We find very little anisotropy between the energies of DAPBs lying in the basal and prism planes. 
Perfect antiphase boundaries in DOig ordered TiaAl are also investigated and their interfacial energies, 
interfacial stresses and local displacements are calculated from first principles through direct supercell 
calculations. Our results are discussed in light of mechanical property measurements and deformation 
microstructure strudies in a Ti-Al alloys. 

1 Introduction 

The presence of chemical short-range order (SRO) in a- Ti-Al was proposed over thirty years ago by Blackburn 
and Williams [1] to explain planar slip exhibited by alloys with concentrations greater than roughly 4 wt. 
% Al (7 at. % Al) [2,3]. Since that time SRO has been discussed as an important factor contributing to 
the anomalous solid-solution strengthening of a-Ti-Al [1, 4-6], as well as the low strain hardening exponents 
exhibited by these alloys [3,4,7]. Only very recently, however, has direct evidence for SRO been provided 
from diffuse neutron scattering measurements [3,8]. These measurements confirm the presence of appreciable 
SRO in Ti-6 wt. %A1 alloys, of a type consistent with the ordered DO19 structure of an TisAl. 

As first proposed by Fisher [9], SRO in solid solutions can provide a significant source of strengthening 
due to the energy cost associated with the creation of "diffuse" antiphase boundaries (DAPB) arising from 
shearing of the lattice by dislocation slip. The DAPB is the analogue for structures with SRO of the 
more familiar antiphase-boundary (APB) in a long-range-ordered crystal structure. In solid solutions with 
chemical SRO, the energy cost of creating a DAPB ("fsROi defined as the excess energy per unit area of 
a DAPB) during slip gives rise to a friction stress on moving dislocations: t sro^ = Isro- A plot of the 
DAPB energy versus the number of shearing events, calculated from first-principles for a-Ti-Al alloys using 
the approach described below, is shown in Fig. 4. This plot is qualitatively similar to theoretical results 
obtained for fec-based alloy systems previously [10-13]. The energy of the DAPB after the first shear (to 
be referred to as 7-^ is seen to be larger than the values (7 2 , 73, etc.) of 75^0 after multiple shearing 
events. As pointed out originally by Cohen and Fine [14], the second slip event in a solid-solution with SRO 
can partially restore the state of order, giving rise to a decrease in Jsro ( c -f- Fig. 4) and a driving force 
for the pairing of leading dislocations in pile-ups. Such pairing has been observed recently by transmission 
electronc microscopy (TEM) for dislocation pile-ups in a- Ti-Al [15]. Fig. 4 also illustrates that after a small 
number of slip events the state of order is largely destroyed and subsequent dislocations experience negligible 
friction stress due to SRO. By this "glide-plane-softening" mechanism, SRO can promote planar slip in solid 
solutions, as proposed for a-Ti-Al alloys by Blackburn and Williams [1]. 



Several TEM-based experimental studies have been undertaken to estimate the magnitude of tsro from 
observed dislocation structures [3, 13, 16-21] in short-range-ordered alloys. To date all of this work has been 
devoted to fee alloys, with the exception of the very recent experimental study by Neeraj et al. [3] for 
hep-based a-Ti-Al. In this study, magnitudes of 7 2 for DAPBs on prism planes were estimated from the 
critical stress required to form screw dislocation dipoles from elongated "hair-pin" loops, while the difference 
7i — 72 was estimated from measured spacings between paired leading edge dislocations in pile-ups. This 
work represents the first direct measurement of the friction stress due to SRO in a-Ti-Al and is an important 
contribution to the detailed understanding of deformation properties in these alloys. The purpose of the 
research reported in this manuscript is to provide an independent theoretical estimate of DAPB energies in 
this system. 

A number of theoretical studies have been undertaken to calculate the magnitude of DAPB energies 
for fee and bec alloys [10-14,22-28]. These studies have employed related methods which require (i) a 
description of the equilibrium state of short-range order in the alloy based upon the SRO parameters, as 
well as (ii) values of the interatomic interaction parameters which parametrize the energy associated with 
changes in atomic configurations. From the SRO parameters the number of favorable and unfavorable bonds 
created by a slip process can be derived and converted to an excess energy per unit area from the values 
of the interaction energies. To date the most comprehensive of such theoretical approaches was developed 
by Schwandcr et al. [12] who formulated general expressions for ~fsRO based upon an arbitrary number 
of SRO parameters. For Ni-Cr alloys, Schwander et al. made use of measured SRO parameters to derive 
the magnitudes of pair interaction energies through the inverse-Monte-Carlo method [29]. This approach 
was employed more recently by Plessing et al. [13] for Cu-Al and Cu-Mn alloys, where good agreement 
was obtained between calculated values and estimates of 7 SRO derived from friction stresses deduced from 
pile-up geometries in post-deformed samples. For a-Ti-Al alloys, measured values of the short-range-order 
parameters are currently unavailable and the theoretical approach of Schwander et al. cannot be applied to 
the calculations DAPB energies in this system at present. An alternative approach relies upon the use of 
quantum-mechanical calculations to derive the values of the interatomic interaction parameters from first- 
principles [30, 31] . These parameters can be employed in statistical-mechanical calculations to compute values 
for short-range-order parameters, allowing estimates of Jsro t° be made without the need for experimental 
input. Such a first-principles approach has been applied to the study of SRO hardening in Cu-Au alloys 
by Mohri et al. [11]. In the present study first-principles calculations of the composition and temperature 
dependencies of Jsro are undertaken for a-Ti-Al alloys. To the best of our knowledge, this work represents 
the first time a calculation of DAPB energies has been undertaken for an hep-based alloy system. 

In the next section we begin by discussing first-principles results for the structure, energetics and stress 
of APBs in ordered a-i T13AI. The APB energy results serve as a basis for testing the predictive power of the 
calculated model interatomic interaction parameters for hep-based Ti-Al alloys. The next section is devoted 
to a detailed discussion of our study of DAPB energies. Our calculated values of 75^0 for DAPBs on prism 
planes are found to be consistent with the recent experimental estimates of Neeraj et al. mentioned above. 

2 Anti-Phase Boundaries 
2.1 Methodology 

The APB energies were calculated using the local density approximation (LDA) to density functional theory 
(DFT), as implemented in the VASP package [32-34], which relies on ultrasoft pseudopotentials [35] and 
a plane wave basis. The plane wave energy cutoff was set to 270 eV while the number of fc-points was 
chosen to be at least 3500/(atom _1 ) for all structures. 1 A supercell geometry was used, where a periodic 
array of APBs is set up. It was verified that the APBs are sufficiently spaced from one another to avoid 
interactions by computing the APB energy for different supercell sizes. Separations of 5 and 16 monolayers 
were considered. Atomic positions were fully relaxed and the dimensions of the supercell perpendicular to 

1 e.g., (3500/n) fc-points would be selected for a structure with n atom per unit cell. 
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the plane of the APB were allowed to relax to give zero normal stress. The lattice parameters along the 
plane of the interface were frozen to the equilibrium lattice parameter of the DO19 structure. 

The error bars we provide were derived by studying the convergence of our results with respect to plane 
wave energy cutoff, number of fc-points, supercell size and the tolerance for residual forces and stresses for 
the numerical determination of the equilibrium geometry. 

2.2 Results and Discussion 

We consider the APB lying in the basal [0001] and prism [llOO] planes, as these are the main active slip 
systems in Ti-Al alloys. The calculated APB energies are reported in Table 1. While there is only one type 
of basal APB, there are two types of prism APB, differing by the location of the APB plane, as illustrated 
in Fig. 1. 

As shown in Table 2, the calculated APB energies are comparable to the result of earlier all-electron 
first-principles calculations [36], indicating that the pseudopotential approximation used provides reliable 
results. Our energies are systematically smaller, however, perhaps due to a more thorough treatment of the 
relaxations made possible by the high computational efficiency of pseudopotential-bascd methods. 

Table 2 also shows that the APB energies calculated from first principles differ markedly from experimen- 
tal estimates [37], which were obtained from measurements of the spacing between dissociated dislocations, 
inferring the interface energy through mechanical equilibrium considerations involving the trade-off between 
clastic and interfacial energy. As discussed in the next section, we estimate that the accuracy of first- 
principles calculations may be responsible for an overestimation of at most 40% of the APB energies. This 
error is significantly smaller that the observed disagreement with experimental estimates. The level of dis- 
crepancy between theory and experiment displayed in Table 2 has been observed for related intcrmctallic 
systems [38,39]. While the source of this discrepancy remains unclear, we note that the equilibrium spac- 
ings implied by our APB energies (using the same elastic constants as in [37]) are 3.0 nm, 7.0 nm, and 1.6 
nm, respectively, for the basal, prism I, and prism II APB. These values are all very close to the minimum 
dislocation spacings measured by [37]. 

The calculated APB energies are consistent with an intuitive bond-breaking picture. The prism I, basal 
and prism II APB respectively create 0,1 and 2 Al-Al bonds, implying that the basal APB should have an 
energy than is approximately equal to the average of the prism I and prism II APB (293 mJ/m 2 ). 2 Our 
calculations also provide interface stress and displacements (see Table 1). An analysis of their sign and 
magnitude reveals interesting features of the relaxation taking place in the vicinity of an APB in TiaAl 
which may help in interpreting the stress and displacement fields that can be quantified with TEM. For the 
two APBs that result in a change of the nearest neighbor environment (basal and prism II), the interface 
contracts perpendicular to the plane of the APB (di < 0) and a tendency to expand along the plane of 
the interface (C22 > and er 33 > 0). This behavior is reversed for the prism I APB, where no change 
in the nearest-neighbor environment occurs. The fact that, in all cases, the stress along the interface and 
the displacement perpendicular to it have opposite signs can be understood from a volume conservation 
argument. However, we were unable to rationalize the relative magnitudes of the stresses and displacements 
on the basis of intuitive atomic size considerations. This is not surprising given that such a picture is unable 
to describe the differences in bond length between hep Ti and (metastable) hep Al. According to our first- 
principles calculations, in these two pure elements, the Ti-Ti bonds are longer than Al-Al bonds lying along 
the (0001) plane while the opposite is true for other nearest neighbor bonds. 

2 The atomic density in the plane of the APB is slightly different in basal and prism APB so that even in an ideal bond-break 
picture, the basal APB energy would differ slightly from the average energy of the prism APBs. 
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3 Diffuse Anti-Phase Boundaries 



3.1 Methodology 

The calculation of diffuse anti-phase boundary (DAPB) energies requires the construction of a thermody- 
namic model of the alloy. The cluster expansion formalism [30,31,40] was used to parametrize the configu- 
rational dependence of the energy of the alloy system: 

E ( (J ) = X/^ i(Ti + X] Jij a i a j + X] JijkO-i&jCrk + ■ ■ ■ (1) 
* {hi} {i,j,k} 

where er, denotes the occupation of site i (<7j = 1 for Al, <7j = —1 for Ti) while the Ji... are effective 
interaction parameters to be determined from first principles and where the summations are taken over 
every singlet, pair, triplet,... of lattice sites. The so-called effective cluster interactions (ECI) Ji... embody 
all the contributions to the energetics of the alloy, including relaxations of the atoms away from their ideal 
lattice sites. These ECI were fitted to the energies of a set of 55 periodic structures calculated from first- 
principles using the same techniques and parameters as for the APB calculations. This procedure was carried 
out using the MAPS software package [41] which automates most of the tasks required to apply the cluster 
expansion formalism. 

The ECI provide all the information we need to compute the DAPB energies. First, the ECI can be 
directly used to perform Monte Carlo simulations [42] that allow a first-principles computation of the phase 
diagram of the alloy system and degree of short-range order present in the alloy at a given temperature and 
composition ([43], and references therein). To this effect, we set up a Monte Carlo simulation cell consisting 
of a 12 x 12 x 6 superccll of the hep unit cell. The initial configuration was chosen to be random and 
15000 equilibration Monte Carlo passes were performed before the short-range order was computed during 
an additional 2500 passes. The size of simulation cell and the number of simulation steps are such that the 
statistical noise in the computed short-range order parameters is orders of magnitude smaller than the error 
arising from the precision of the ECI themselves. 

The DAPB energies Edapb per unit area can then be directly expressed in terms of the ECI and the 
short-range order parameters obtained from the Monte Carlo simulation: 

Edapb = ^ ^ (( a i a j) ~ ( a s(i) a s(j))) J ij + ■■■ ( 2 ) 
{id} 

where S is a function that maps a given site i onto its new position, after the anti-phase boundary has been 
created and where (<Ji<Jj) denotes the thermodynamic and spatial average of the correlation between the 
occupation of sites i and j. The sum is taken over all pairs (triplets, etc.) that are not equivalent under 
a translation lying in the DAPB plane which leaves the lattice invariant. A is the area of a repeat unit 
of the plane of interest. This formulation, which is equivalent to the one introduced by Schwander [12], 
has the advantage that it makes the most efficient use of all the information provided by the Monte Carlo 
simulations. It allows us to obtain DAPB energies without having to average over various DAPB energy 
obtained by running a dislocation through various snapshots of the Monte Carlo cell during the simulation. 
This averaging is implicitly taken into account through the averaging of the correlations (<Ji<7j) over the 
Monte Carlo steps and over all symmetrically equivalent clusters. As the cluster expansion does not directly 
provide information regarding the area A, we assume the disordered phase to be an hep phase with the idea 
c/a ratio and with an average atomic volume equal to the one of the DO19 structure. For a system with a 
small size mismatch (3%, in our case), the error in A introduced by this approximation will amount to at 
most 6%. 

Since we have access to the energy of a finite number of structures, we can only determine a finite 
number of ECI and the expansion defined by Equation (1) needs to be truncated. The choice of the optimal 
number of terms to keep was made through the following statistical analysis. As our goal is ultimately to 
predict DAPB energies, we seek to construct a cluster expansion that exhibits a good predictive power for 
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alloy configurations similar to perfect APB. This exercise docs not merely consists in trying many different 
ranges of interaction until one happens to predict APB energies that agree with the calculated ones. The 
agreement could be fortuitous and the predictive power of the cluster expansion for configurations that arc 
slightly different from a perfect APB, such as a DAPB, could be poor. Instead, we calculate the statistical 
variance of the predicted APB energies for various choices of range of interaction, and look for the choice 
that minimizes that variance. The statistical variance is computed as follows. Let J denote the vector of 
all the selected ECI and let V be their estimated covariance matrix. 3 The energy of a given APB, Eapb is 
related to the vector J of all ECI by Equation (2), where the correlations {piOj) are calculated for a perfect 
DO19 structure. This linear relationship can be written as, Eapb = C T J, where C is some vector. As a 
result, the variance matrix of the predicted APB energy is then simply given by: C T VC. 

3.2 Assessment of the Cluster Expansion 

The ECI of the cluster expansion that were found to minimize the sum of the statistical variances of the 
predicted APB energies are listed in Table 3. The ECI we calculated are very close to the ones that were 
determined in [44] by a fit to the experimental Ti-Ti 3 Al phase boundary. The resulting cluster expansion 
predicts APB energies that agree remarkably well with the calculated APB (see Table 4), given that the 
perfect APB energies were not used in the fit of the cluster expansion. The difference between the calculated 
and predicted energies is of the same order of magnitude as the error bars, confirming that our variance 
calculations indeed provide a very good measure of the true predictive power of the cluster expansion. A 
similar variance analysis will be used in the next section to estimate the errors on our calculated DAPB 
energies. 

Although the selected range of interaction is by construction the one that will best predict the APB 
energies, it is still possible that longer range interactions may have an important effect on the calculated 
thermodynamic properties of the alloy. To investigate this concern, longer range interactions were gradually 
added. Adding up to 8 additional pair interactions and 6 triplet interactions changed the temperature of the 
DO19 ~^hcp transition by less than 200K and we are thus confident that our short-range cluster expansion 
is sufficiently accurate for our purposes. 

Table 5 gives a few characteristics of our cluster expansion. The precision of the cluster expansion is 
quite good, as can be seen by the value of both the mean square error and the cross-validation score. 4 

The temperature T c of the DO19 -^hep transition was calculated to be 2000K, which is about 40% 
higher than the known experimental transition temperature of 1420K. Such an ovcrestimation of transition 
temperatures is quite common in first-principles calculations based on the local-density approximation and 
short-ranged cluster expansions (a list of systems where such first-principles phase diagram calculations 
have been undertaken can be found in [30]). To correct for this overestimation of T c , we normalize all our 
temperatures by T c . Experimental measurements should then be compared with the calculations performed 
at the same normalized temperature. In Fig. 2, the Ti-rich portion of the calculated phase diagram is seen 
to agree very well with the experimentally measured one, when both are normalized to have the same T c . 

The error in the calculated transition temperature provides us with an estimate of the error introduced 
by (i) the use of the Local Density Approximation in the first-principles calculations (ii) neglecting non- 
configurational sources of entropy, such as lattice vibrations and electronic entropy, and (iii) inaccuracies 
in the description of the configurational free energy arising from the truncation of the cluster expansion to 
a finite number of terms (see above). We thus expect our perfect and diffuse APB energies to be biased 
upward by at most 40% (this correction has neither been applied to our result nor been included in the error 
bars, in order to allow comparison with other first-principles results). 

3 Let Y be an n-dimcnstional vector and X be a n X k matrix. From standard statistical analysis, the covariance matrix of a 
parameter vector b obtained from a least square fit mirif, {(Y - Xb) T (Y - can be estimated by (X T X) 1 a 2 where a 2 
is the mean square error of the fit. 

4 The cross-validation score is defined to be I ^™ =1 I _B S — E s J I , where E s is the calculated energy of structure s 
and E s is the value of the same energy predicted using the cluster expansion fitted to the remaining n — 1 structures. 
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3.3 Results and Discussion 



Our first main result is that Ti-rich Ti-Al alloys do indeed exhibit a substantial amount of short-range 
order, as shown in Fig. 3, which plots the Warren-Cowley short-range order parameters [45] as a function 
of distance. This short-range order is of the DOig type. Figures 4 and 5 show the calculated values of the 
DAPB energies r y k for the prism DAPB for the temperatures and concentrations marked by a cross in Fig. 
2. Energies of the basal DAPB are not shown, as they are nearly identical to those for the prism DAPB, as 
illustrated in Fig. 6 for one specific temperature and concentration. In these plots, the abscissa k represents 
the magnitude of the shift (in multiples of lattice vectors) between the two half crystals on each side of the 
DAPB, or, alternatively, the number of dislocations needed to create the DAPB. Error bars are not shown 
for clarity, but amount to about 15% of the calculated DAPB energies resulting from the statistical noise in 
the fitted ECI. 

Note that the two types of prism APB can no longer be distinguished in a short-range ordered alloy 
because the DAPB plane cuts through various short-range ordered domains which have a random relative 
phase, resulting in a single average interfacial energy. The plots exhibit the characteristic oscillations that 
are expected from the observation that, in a perfectly ordered DOig structure, every other dislocation would 
restore a perfect state of order. In a short-range ordered alloy, the ability of subsequent dislocations to 
restore the original state of order decreases with the number of dislocations as a result of the short-range 
nature of the state of order. 

Our estimates of the DAPB energies are comparable to those found experimentally [3]. For a Ti-6 wt. 
% Al alloys annealed at temperatures between 600 and 350°C, experimental estimates of 7 2 ranged between 
25 and 48 mJ/m 2 , while a value of 6 mJ/m 2 was quoted for the difference 7i — 7 2 - Our results indicate that 
for an alloy of the same composition (10 at. % Al) and for the same range of relative temperature (T/T c 
ranging from 0.45 to 0.6), j 2 = 15 mJ/m 2 while r y 1 — 7 2 ranges from 7 mJ/m 2 to 10 mJ/m 2 . These values 
of 7 2 arc within a factor of 2 relative to experimental results, while our estimates of the difference j x — 7 2 
are comparable. 

An important finding is that for a substantial portion of the region of solid solubility of Al in Ti, the 
DAPB energies r y 1 are larger than 2 — 4 mJ/m 2 , values which have been quoted [46] as the thresholds above 
which transitions from homogenous to planar slip occur in fcc-based Cu-Zn, Cu-Al, and Cu-Mn alloys. This 
result supports the proposition that the presence of SRO contributes to the predominantly planar slip in 
Ti-Al [1]. 

The fact that prism and basal DAPB have an almost identical energy is again consistent with a simple 
bond-breaking picture. 5 Despite the small anisotropy in the average DAPB energies, the behavior of each 
type of DAPB under stress may differ. The local DAPB energy density may very well fluctuate as wc 
move along the interface. For instance, a prism DAPB may contain portions that are locally analogous to 
a low energy prism I APB and portions that are locally similar to a high energy prism II APB, while such 
fluctuations should be much smaller in basal DAPB where only one type of local order exists. The average 
DAPB energy we compute is the quantity of interest in order to determine, for instance, the total work 
required to deform the material by a given amount. However, for properties such as the critical resolved 
shear stress (CRSS), the average may not provide enough information. A more detailed investigation of the 
spatial fluctuations in DAPB energy density would be needed, a topic we leave for future research. 

Our calculated "f k are of the same order of magnitude as earlier theoretical studies which focused on fee 
alloys [10-13]. One notable qualitative distinction is the fact that, in hep Ti-Al alloys, the second dislocation 
restores the original state of short-range order to a greater extent than in these fee alloys, as can be seen 
by the larger drop between ■y 1 and 7 2 . We attribute this difference to the fact that, in all the fee alloys 
studied, the thermodynamically stable ordered phase in the vicinity of the short-range ordered region is a 
DO22 structure. Due to the tetragonal symmetry of this phase, different short-range ordered domains can 
have distinct orientations. Since the DAPB cuts through these different domains and since the DO22 has 
different periodicity in different directions, it is impossible for the second dislocation to restore the local 

5 This is not an artifact of the use of a short-range cluster expansion. The inclusion of second nearest-neighbor interactions 
allows the basal and prism DAPB energies to differ. 
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order in all domains simultaneously. In contrast, in the case of Ti-Al, the hep-based DO19 structure has the 
same periodicity along any [1100] direction and it follows that the second dislocation can partially restore 
the original local order in all domains simultaneously. 

4 Conclusion 

Our first main result is the confirmation, on the basis of first-principles calculations, that DOi 9 -type short- 
range is present in the hep a-Ti phase of the Ti-Al alloy. We have also calculated the energy of diffuse anti- 
phase boundaries (DAPB). We find that for an alloy at a composition of 10 at. % of Al and a temperature 
between 600 and 300°C, a single dislocation creates a DAPB with an interfacial energy ranging from 22 
to 25 mJ/m 2 while a second dislocation creates a DAPB having an interfacial energy of 15 mJ/m 2 . These 
values are significantly higher than 4 mJ/m 2 , a value estimated to be sufficient to induce a transition from 
homogeneous to planar slip in fec-based alloys [46]. 

The computer programs used to complete this work can be obtained by contacting the authors. These 
include an easy-to-use Monte Carlo code that provides, among other quantities, short-range order parameters 
and automatic phase boundary plotting capabilities. A companion code implements the Schwander et al. 
[12] formalism to its full generality, allowing the calculation of DAPB of any orientation, for any lattice type, 
and for any number of short-range order parameter. 
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Tables 





basal 


prism I 


prism II 


7 (mj/m 2 ) 


249 ±5 


108 ±2 


477 ±5 


er 22 (J/m 2 ) 


1.4 ±0.2 


-1.5 ±0.2 


3.8 ±0.2 


CT33 (J/m 2 ) 


0.8 ±0.2 


-0.1 ±0.2 


0.3 ±0.2 


di (A) 


-0.02 ±0.01 


0.16 ±0.01 


-0.10 ±0.01 


d 3 (A) 


-0.15 ±0.03 









Table 1: Calculated Interface Energy 7, Interface Stress a and Displacement d of the Basal and Prism APB. 
Stresses and displacements are reported in cartesian coordinate system where the two first basis vectors are 
respectively defined by the unit normal and the Burgers vector of the APB. Compressive stresses are positive. 
The negative sign of d% indicates that Al-Al bonds at the interface are shortened by the displacements. All 
unreported stesses and displacements are zero by symmetry. 





7 (mJ/m 2 ) 




basal 


prism I 


prism II 


present 


249 


108 


477 


first-principles [36] 


300 


133 


506 


experiment [37] 


63 


42 


84 



Table 2: Comparison between calculated and Measured Interface Energy 7 of the Basal and Prism APB. 



Pair 


Interaction (eV) 


Rcf [44] 


(0,0,0)- a, £) 


0.057 ±0.005 


0.0588 


(0,0,0) -(1,0,0) 


0.047 ±0.003 


0.0588 


(0,0,0)- (l-hi) 


-0.014 ±0.004 


-0.0235 



Table 3: ECI of the Cluster Expansion. Coordinates are expressed in the conventional hexagonal coordinate 
system. Error bars are one standard deviation wide. 
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7 (mJ/m 2 ) 


basal 


prism I 


prism II 


supercell 


249 ±5 


108 ±2 


477 ±5 


clus. exp. 


234 ± 34 


71 ±20 


427 ±55 



Table 4: APB Energies Obtained from Supercell Calculations and Predicted from a Cluster Expansion. 
Error bars on the predicted energies are one standard deviation wide. 



Number of structures: 


55 


Atom/unit cell: 


2-12 


Mean square error: 


0.031 eV/atom 


Cross-validation score: 


0.034 cV/atom 


Predicted T c of DOi 9 -^hcp: 


2000K±25K 


Experimental T c of DO19 -^hcp: 


1420K 



Table 5: Characteristics of the Cluster Expansion. 
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Figure 1: DO19 structure and the location of prism I (dashed line) and prism II (solid line) planes. The size 
of the circles distinguishes distinct (0001) planes of atom. Al atoms are shown in white, Ti atoms in black. 
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Figure 2: Calculated and Experimental Phase Boundaries (from [47]). Also shown are the points where the 
DAPB energies were calculated. 
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SRO parameter at x A1 =10% and T/T c =0.60 
0.15 n < ' ' ' 1 




Figure 3: Calculated Warren-Cowley short-range order parameter as a function of distance (normalized by 
nearest neighbor distance). 



y k for prism DAPB at x A1 =12% 




Figure 4: Calculated Prism DAPB Energies as a Function of Temperature. 



y k for prism DAPB at 17T C =0.60 
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Figure 5: Calculated Prism DAPB Energies as a Function of Concentration. 
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Figure 6: Calculated Prism and basal DAPB Energies. 
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Figure Captions 

Figure 1. DOi 9 structure and the location of prism I (dashed line) and prism II (solid line) planes. The size 
of the circles distinguishes distinct (0001) planes of atom. Al atoms are shown in white, Ti atoms in black. 

Figure 2. Calculated and Experimental Phase Boundaries (from [47]). Also shown are the points where the 
DAPB energies were calculated. 

Figure 3. Calculated Warren-Cowley short-range order parameter as a function of distance (normalized by 
nearest neighbor distance). 

Figure 4. Calculated Prism DAPB Energies as a Function of Temperature. 

Figure 5. Calculated Prism DAPB Energies as a Function of Concentration. 

Figure 6. Calculated Prism and basal DAPB Energies. 
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